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ABSTRACT 

It is well known that the mere location of a Beat Cepheid model in a Period Ratio vs. Period diagram (Petersen 
diagram) puts constraints on its metallicity Z. But these bounds are sensitive to the mixture of elements that 
are lumped into the parameter Z. In this short paper we update the previous results that were based on the 
Grevesse-Noels solar mixture to the recent, revised Asplund, Grevesse & Sauval (2005) (AGS) solar mixture. 
We also examine the effect of the envelope depth on the accuracy of the computed pulsation periods. We find 
that for low period Cepheids with high Z the customary approximation of envelope pulsation breaks down. 
It is necessary to compute stellar models that extend to the center and to include burning and composition 
inhomogeneities in the modeling. Fortunately, however, most Beat Cepheids that have been observed so far 
seem to avoid that regime. 

Subject headings: (stars: variables:) Cepheids, stars: oscillations, stars: rotation, galaxies: abundances 



1. INTRODUCTION 

A recent paper by Buchler & Szabo (2007) (Paper I, here- 
after) described the results of a modeling survey of Beat 
Cepheids that pulsate in the fundamental (F) mode and the 
first overtone (Ol) simultaneously. In order to allow a con- 
venient comparison with the observations, Paper I presented 
the results in a Pio vs. Pq diagram, a so called Petersen di- 
agram, as a function of metallicity Z, where Pq and Pi are 
the fundamental and first overtone periods, respectively, and 
Pio = Px/Pq, The two periods can easily be extracted from 
the observed light curves. From the metallicity dependent Pe- 
tersen diagram one can set bounds on the metallicity Z of the 
star. This way of determining the metallicity is independent of 
and complementary to all other stellar and galactic metallicity 
probes. 

Buchler & Szabo (2007) found that the Petersen diagrams 
are rather insensitive to stellar rotation rates and to the galac- 
tic helium enrichment function Y = Y(Z) that they used. 
However, the results exhibited sensitivity to the chemical mix- 
ture that is lumped into the metallicity parameter Z, mainly 
through their effect on the opacities. The survey of paper I 
used the OPAL GN93 opacities (Iglesias & Rogers (1996) for 
the solar mixture of Grevesse & Noels (1993), merged with 
the Alexander & Ferguson (1994) opacities at the lower tem- 
peratures (log T <3.95). 

In the last couple of years Asplund, Grevesse & Sauval 
(2005) have dramatically revised the solar chemical compo- 
sition of metals. The most important change is a decrease of 
elements with atomic numbers Z < 11, but dominated by O, 
in favor of elements with Z > 1 1 of which Fe has the most 
important effect on the opacities. In view of the sensitivity 
of the Petersen diagrams to the chemical mixture we have 
thought it useful to update our survey with OPAL opacities 
that reflect these new abundances. 

In this study, as in Paper I, we use a solar elemental mixture 
even in low metallicity Cepheids out of expediency, because 
there is no observationally derived information available that 
would let one improve much upon such an approximation at 
this time. 
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Fig. 1. — Helium galactic enrichment function Y(Z) that is used in 
the model computations. 

In a linear study such as this the Beat Cepheids are defined 
as being simultaneously unstable in the fundamental (F) and 
the first overtone (Ol) modes. However, it is well known from 
nonlinear Cepheid studies {e.g., Kollath et al. (2002); Kollath, 
Beaulieu, Buchler, Yecko (1998)) that not all stars in this 'lin- 
ear Beat Cepheid instability strip' will actually undergo sta- 
ble Beat pulsations (a large fraction will be either F or Ol 
Cepheids). This computationally convenient, linear definition 
therefore considerably overestimates the 'Beat Cepheid insta- 
bility strip' (broader by at least a factor of 3 - 4 in logX^ 
in a theoretical HR diagram (Szabo & Buchler 2008). This 
in turn increases considerably the separation of the upper and 
lower boundaries in the Petersen diagrams with a concomitant 
increase in the uncertainty of the stellar metallicities. The lin- 
ear properties of stars do not allow one to infer the modal 
selection which is intrinsically nonlinear (Stellingwerf 1974; 
Buchler & Kovacs 1986). Therefore, for a tighter, and correct, 
definition and study of Beat Cepheids it is necessary to resort 
to a full amplitude hydrodynamic survey of models which re- 
quires a huge computing effort. Such a survey will be pre- 
sented in a separate paper (Szabo & Buchler (2008)). 

In §2 we display the computed Petersen diagrams for the 
AGS solar elemental composition. In §3 we discuss the ad- 
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Fig. 2. — Pi / Po VS. logPo (Po on top axis) where Po and Pi are the fundamental and first overtone periods. The lines delimit the ranges for 
which both F and 01 are linearly unstable: Lower boundaries in the left, and upper boundaries in the right panels. The metallicity increases 
downward in the figure from Z = 0.001 (top line) to 0.016 in steps of 0.001, and then from in steps of 0.002 to 0.020 (bottom line). (In the 
left panels the lowest Z values are not shown to avoid overcrowding). The long-dashed line refers to the specific AGS composition of the Sun, 
Zq — 0.0122 and Yq = 0.2486. For reference we have superposed the location of most known Beat Cepheids. 
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TABLE 1 

Galactic Metallicities derived from the Beat 
Cepheids for the AGS and GN93 elemental mixtures. 

Z m i n and Zmax represent the lowest and largest Z values that are obtained 
for all Beat Cepheids in the respective galaxy. (Z) is the average of the 
centroids of the mininum and the maximum Z of each star in the galaxy. 
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equacy of envelope models for Cepheids, and find that it is 
necessary to compute full stellar models for the short period, 
low metallicity Cepheids. We conclude in §4. 

2. PETERSEN DIAGRAMS WITH THE AGS CHEMICAL 
COMPOSITION 

This survey proceeds in parallel with the one presented by 
Buchler & Szabo (2007) for the GN93 solar chemical mixture, 
and we do not repeat a description of the procedure and of the 
modeling details here. We just mention that we have kept 
the same helium enrichment function Y(Z) (Fig 1) that was 
based on an interpolation between the 'standard' values for 
the Galaxy and the two Magellanic Clouds. We stress again 
that the results are very broadly insensitive to this relation. 
We have used the OPAL website to generate opacities as a 
function of metallicity Z with the mixture of elements that 
are taken from the AGS table of solar abundances. 

Our newly computed Petersen diagrams for the AGS mix- 
ture are displayed in the top panels of Fig. 2. The left-side 
and right-side panels show the lower and upper bounds, re- 
spectively, of the allowed Beat Cepheid regions as a function 
of metallicity Z. We alternate the line-type (or color) for the 
successive Z values to make the diagrams easier to use. In 
order to avoid crowding we have omitted the labels 0.001 to 
0.005 on the upper curves in the left-side panels. 

The dashed curves in the top panels in Fig. 2 show the re- 
sults for the opacity with the specific AGS composition for the 
Sun, viz. Z Q = 0.0122 and Y Q = 0.2486. Despite the large 
difference in Y (this point lies below the range of the diagram 
in Fig. 1), this curve falls close to the 0.012 lines, illustrating 
again the insensitivity to the exact values of Y. 

For an easy comparison we have also juxtaposed the GN93 
composite Petersen diagram of Paper I in the bottom panels 
of Fig. 2. The GN93 models have been recomputed with 
the same zoning (200 mass-shells) as for the new AGS re- 
sults, which is a little finer than that used in Buchler & Szabo 
(2007). It is noteworthy that the AGS mixture causes an over- 
all downward shift of the curves, especially for the shorter 
periods, and a downward stretching for the longer periods. 
One therefore obtains generally lower metallicities with the 
AGS solar mixture than with the GN93 solar mixture. This is 
of course expected since AGS have revised the metallicity for 
the Sun down to Z G) = 0.0122. 
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Fig. 3. — The period ratio Pw as a function of the temperature of 
the envelope depth, T c (in 10 6 K), for the sequences of models of 
Table 1 with the indicated fundamental periods (in days). 



The locations of the known Beat Cepheids are superposed 
in Fig. 2 as in Buchler & Szabo (2007), with open squares 
for the Galaxy, x's for M33, filled triangles for the LMC and 
crosses for the SMC. We have limited the figure to periods 
longer than one day despite the fact that a few such stars are 
known in the SMC. The reason will be given in §3. 

The diagrams of Fig. 2 are used as follows: If a Beat 
Cepheid falls between the upper and lower boundaries for a 
given Z value, that means that models with that value of Z 
are compatible with the observed period constraints for the 
given Beat Cepheid. By finding all the compatible Z values 
one can set upper and lower bounds on Z. In practice, one in- 
terpolates between the displayed curves in order to refine the 
Z bounds. 

In Table 1 we present the metallicities for the 4 galaxies 
in which a significant number of Beat Cepheids are known. 
From the composite Petersen diagram we infer for each star 
(i) the range of possible Z, i.e. Z m i n [ and Z max ,i which, as we 
have seen, represent the uncertainty range of Z for each Beat 
Cepheid (due to the width of the instability strip). We define 
■^min as the minimum over the galactic sample of these Z m i n ,i 
and Z max as the maximum of the Z maX:1 . Finally we introduce 
the galactic average (Z) as the average of the centroids of the 
ranges \ {Z min ^ + Z max ,i). 

With the new opacities we obtain an average metallicity of 
(Z) = 0.0088 and a range of uncertainty of 0.0055 to 0.0133 
for the Galactic Beat Cepheids. This is to be compared with 
an average (Z) = 0.01 18 and a range of 0.0073 to 0.0182 with 
the GN93 opacities. As noted in Paper I the GN93 average 
(Z) from the Galactic Beat Cepheids is quite low compared 
to the GN93 value for the Sun Z Q = 0.017, the but the AGS 
the average (Z) of the Beat Cepheids is now much closer to 
the AGS value of Z Q =0.0122. 

We note again here that the large, indicated spread Z m - m 
to Z max represents a combination of the actual spread in the 
galactic metallicity, and of the uncertainty that arises from the 
use of the linear definition of Beat Cepheids. As mentioned 
earlier, we hope to reduce the latter with a survey of full am- 
plitude Beat Cepheid models Szabo & Buchler (2008). 




TABLE 2 

Model Sequences used in the Tests of §3 
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Fig. 4. — (a) Model structure, logp[g/cm 3 ] vs. logT[K] for the 
5 models of Table 2 on the fundamental blue edge with fundamental 
periods, 0.5, 1.0, 2.0, 4.0 and 7.0 d; (b) Inverse sound speed (in units 
of 10 cm/s) vs. radius [cm]; (c) Inverse sound speed vs. relative 
radius. 



3. DISCUSSION OF THE MODELING UNCERTAINTIES 
- INADEQUACY OF ENVELOPE MODELS 

Ever since the early days of Cepheid modeling, it has been 
customary to approximate Cepheids with envelope models of 
uniform composition, because the pulsation is largely con- 
fined to the outer parts of the stars. One thus imposes a rigid 
core of radius R c at some temperature T c , located outside the 
burning region where the luminosity L c is constant. This ap- 
proximation has been justified on the one hand, because the 
radial displacement eigenvectors for the low lying pulsation 
modes decay rapidly inward and one therefore sets Sr = at 
r = R c , and, on the other hand, because the computed peri- 
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ods are found not to change much when one varies the depth 
of the envelope models. However, for the Petersen diagrams 
we require a higher relative accuracy (< 0.001) in the periods, 
as we now discuss. 

A reexamination of this assumption shows a breakdown for 
models with low periods and high Z. It turns out that in this 
regime the periods, and concomitantly the period ratios, ex- 
hibit an uncomfortable dependence on R c or T c . Furthermore, 
for the lowest periods there is a poor leveling off even with 
T c up to (perhaps unreasonably large) values of 5.5xl0 6 K 
which is beyond the temperature up to which one can neglect 
nuclear burning and spatial composition inhomogeneities. 

This T c dependence is illustrated in Fig. 3, which plots 
the period ratios Piq vs. T c for a sequence of 5 models 
(see Table 2) that have been chosen with a composition of 
X = 0.7045 and Z = 0.016 and that are all located at the 
high temperature side of the region where both the fundamen- 
tal and the first overtone modes are linearly unstable. The 
fundamental periods of the sequences are Po = 0.5, 1.0, 2.0, 
4.0, and 7.0d, respectively. The shape for the leveling off is 
found to be the same for all sequences, but while the variation 
in Pio remains within an acceptable range for the long period 
models, this is not the case for the short period models. For 
example, in the 4 d sequence a change from T c = 2.5 x 10 6 to 
3.0 x 10 6 K leads to an increase of P w of less than 0.0002, but 
for the 0.5 d sequence the same change results in an increase 
of Pio of 0.0035. This is a large amount if one considers the 
Petersen diagram of Fig. 2. Furthermore, for the low period 
models the leveling off is not complete even by the time one 
reaches T c = 5.5 x 10 6 K. 

In Fig. 4 we present some other relevant properties of the 
sequence of models from Table 2. Thus the top panel displays 
the models in the log p vs. log T plane, and reveals a very sim- 
ilar structure for the 5 models and demonstrates that nothing 
special is occurring for the 0.5 d models. 

The middle (b) and bottom (c) panels exhibit the inverse of 
the sound speed vs. radius, r, and vs. relative radius (r/P*), 
respectively, where P* is the equilibrium stellar radius. The 
dots correspond to the location of the rigid inner core P c , de- 
fined to lie at T c — 2.5 x 10 6 K in these models. In actual 
radius, the low period models are much deeper, but this is not 
the case in relative radius. The reason for presenting these 
quantities is that the fundamental period Po is roughly pro- 
portional to the stellar sound traversal time 
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The combined effect of the high sound speed in the stellar 
core and the decay to zero of the eigenvectors cause the con- 
tribution from the neglected penetration of the core to be very 
small in general. It is clear from Fig. 4(b) that the contri- 
bution to the period from a deepening of the envelope by 




Fig . 5 . — Sensitivity to the depth of the envelope T c of the Pi /Pa vs. logPo [d] Petersen diagrams (for an AGS mixture). The solid lines show 
the results for an envelope of depth T c = 2.5 x 10 6 K and the dashed curves for T c = 4.5 x 10 6 K. Same notation as in Fig. 2. 



AR C , namely AR C / (C s ), is larger for the long period mod- 
els. However, when one considers it relative to the period, 
one finds a much larger contribution for the envelope mod- 
els with the smaller periods. For example, from the figure 
we infer that in the 7.0 d model a small change of 10 10 cm in 
R c , i.e. AR C /R* ~ 0.003, would cause an increase of Po by 
~ 0.005 d, i.e. a ~ 0.07 % relative period change. In the 0.5 d 
model the same change of R c , i.e. AR C /R* ~ 0.02, would 
cause a considerable ~1 % change in the period. 

Fortunately, an accuracy in the pulsation periods of a per- 
cent or less is generally obtainable and sufficient for most pur- 
poses, so that envelope models are therefore quite adequate. 
Here, however, one puts higher demands on the Cepheid mod- 
els, because one wants to extract a Cepheid's metallicity from 
its location in the Petersen diagrams, for which a relative ac- 
curacy of ~ 10~ 3 is necessary. 

In Fig. 5 we present a composite Petersen diagram that 
compares Beat Cepheid models in which the depth is set at 
our standard T c = 2.5 x 10 6 K (solid lines) with those at 
T c = 4.5 x 10 6 K (dashed lines). One notices that the 
most unreliable models are to be found in the high Z, low 
period domain. Luckily, the vast majority of the observed 
Beat Cepheids avoid the regime in which the uncertainties are 
largest. 

Fig. 5 shows that one would be ill advised to use envelope 



models for the 1.0 d stars that are found in the SMC (Mar- 
quette et al. (2007)). For these Beat Cepheids it is necessary to 
use models that extend deeper, if not full stellar models. This 
is, of course, a much harder problem, in that it requires the in- 
clusion of nuclear burning and of the spatial composition vari- 
ations that occur during stellar evolution. In other words, one 
needs a stellar evolution code that incorporates a linearization 
to compute the fully nonadiabatic pulsation frequencies with a 
time-dependent mixing length recipe to capture the pulsation 
- convection interaction. 

4. CONCLUSION 

We have computed the Beat Cepheid models that make up 
the composite Petersen diagrams as a function of metallicity 
Z with the Asplund, Grevesse & Sauval (2005) solar mixture. 
This change causes a noticeable reduction of ~ 25 - 30 % in 
the average (Z) of the four galaxies in which a sizable num- 
ber of Beat Cepheids are known. It considerably improves 
the agreement between the average Galactic metallicity and 
the Sun's metallicity over what was obtained with the GN93 
abundances. 

While the elemental mixture of the low Z stars is most 
likely not quite solar, this is the best approximation we can use 
in the absence of detailed determinations of chemical compo- 
sitions from high resolutoin spectra of such stars. 
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Because we needed to construct Beat Cepheid models down 
to periods as short as 0.5 d as have been found in the SMC 
by Marquette et al. (2007), we have been led to reexamine 
the accuracy of commonly used envelope models for Cepheid 
stars. We find that because of the high accuracy demands 
of Petersen diagrams such envelope models become inade- 
quate for short period, high Z Cepheids. Good fortune has 
it though that the higher Z Beat Cepheids in the Galaxy and 
in M33 have longer periods, and that the Beat Cepheids with 
smaller periods are observed in the lower metallicity Mag- 
ellanic Clouds. Thus the vast majority of the known Beat 
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Cepheids avoid to some extent the troublesome region of un- 
certainty of envelope models. 
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